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Abstract 

A one-dimensional quantum version of the frustrated XY (planar rotor) 
model is considered which can be physically realized as a ladder of Josephson- 
junctions at half a flux quantum per plaquette. This system undergoes a 
superconductor to insulator transition at zero temperature as a function of 
charging energy. The critical behavior is studied using a Monte Carlo transfer 
matrix applied to the path-integral representation of the model and a finite- 
size-scaling analysis of data on small system sizes. Depending on the ratio 
between the interchain and intrachain couplings the system can have single 
or double transitions which is consistent with the prediction that its criti- 
cal behavior should be described by the two-dimensional classical XY-Ising 
model. 



1 



The two-dimensional frustrated classical XY model has attracted considerable atten- 
tion recently0i'ii'illl'i. It can be related to Josephson-junction arrays in a magnetic field, 
where it is expected to describe the finite-temperature superconductor to normal transition 
in arrays with half a flux quantum per plaquettei. At low temperatures where capacitive 
effects dominate, the array undergoes a superconductor to insulator transition as a function 
of charging energy0'0SEl. These charging effects arise from the small capacitance of the 
grains making up the array and leads to strong quantum fluctuations of the phase of the su- 
perconducting order parameter. The critical behavior is now described by a two-dimensional 
frustrated quantum XY model with a Hamiltonianlll 

H = ~ f E ijf) - E E„' COs(6 r - Or,) (1) 
r \ r / <rr'> 

The first term in Eq. (1) describes quantum fluctuations induced by the charging energy 
E c = 4e 2 /C of a non-neutral superconducting grain located at site r, where e is the elec- 
tronic charge and C is the effective capacitance of the grain. The second term is the usual 
Josephson-junction coupling between nearest-neighbor grains. 9 r represents the phase of 
the superconducting order parameter and the couplings E rr i satisfy the Villain's 'odd rule' 
in which the number of negative bonds in an elementary cell is oddB In a square lattice 
this can be satisfied, for example, by ferromagnetic horizontal rows and alternating ferro- 
magnetic and antiferromagnetic columns of bonds. This rule is a direct consequence of the 
constraint that, for the half flux case, the line integral of the vector potential due to the 
applied magnetic field should be equal to n in units of the flux quantum. 

In this work we consider a one-dimensional frustrated quantum XY (ID FQXY) model0 
which can be regarded as the simplest ID version of the model (1) consisting just of a single 
column of frustrated plaquettes as indicated in Fig. 1. This can be physically realized 
as a periodic Josephson-junction ladder at half a flux quantum per plaquetteffl. In the 
classical limit (E c = 0), the ground state of the ID FQXY model has a discrete Z2 symmetry 
associated with an antiferromagnetic pattern of plaquette chiralities Xp = ilj as indicated 
in Fig. 1, measuring the two opposite directions of the super-current circulating in each 



plaquette. For small E c , there is a gap for creation of kinks in the antiferromagnetic pattern 
of Xpi an d the ground state has long-range chiral order. 

Within a path-integral approach^ , it can be shown that the effective action describing 
quantum fluctuations in the ID FQXY model leads to two coupled XY models in two (space- 
time) dimensions which is expected to have a critical behavior in the universality class of 
the 2D XY-Ising model@ defined by the classical Hamiltonian 

[3H = - J2 iM 1 + a rVr' COs(0 r -6 r ')+ Co r O r ,\ (2) 

<rr'> 

The phase diagram of this model consists of three branches, in the ferromagnetic region. 
One of them corresponds to single transitions with simultaneous loss of XY and Ising order. 
Further away from the branch point, this line of single transitions becomes first order. The 
other two lines correspond to separate XY and Ising transitions. The ID FQXY model 
is represented by a particular path through this phase diagram which will depend on the 
ratio E x /E y between the interchain and intrachain couplings of the model. In this work we 
describe the results of a finite-size scaling analysis of extensive calculations on the ID FQXY 
mo which shows that, in fact, the XY and Ising-like excitations of the quantum model 

decouple for large interchain couplings, giving rise to pure Ising model critical behavior for 
the chirality order parameter. As a consequence, the universality class of the superconductor- 
insulator transition in the related Josephson-junction ladder should then depend on the ratio 
between interchain and intrachain Josephson couplings. 

To study the critical behavior of the ID FQXY model, we find it convenient to use 
an imaginary-time path-integral formulation of the modefi In this formulation, the one- 
dimensional quantum problem maps into a 2D classical statistical mechanics problem where 
the ground state energy of the quantum model of finite size L corresponds to the reduced 
free energy per unit length of the classical model defined on an infinite strip of width L along 
the imaginary time direction, where the time axis r is discretized in slices At. After scaling 
the time slices appropriately in order to get a space-time isotropic model one obtains a 
classical partition function where the parameter a = (Ey/Ec) 1 ^ 2 plays the role of an inverse 
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temperature in the 2D classical model. The scaling behavior of the energy gap for kink 
excitations (chiral domain walls) of the ID FQXY corresponds to the interface free energy 
of an infinite strip in this classical model. For large a (small charging energy E c ), there 
is a gap for creation of kinks in the antiferromagnetic pattern of Xp an d the ground state 
has long-range chiral order. At some critical value of a, chiral order is destroyed by kink 
excitations, with an energy gap vanishing as \at — a c \ u ', which defines the correlation length 
exponent v. Right at this critical point, the correlation function decays as a power law 
< XpXp' > = \p ~ P'\~ v with a critical exponent rj. The free energy per unit length f(a) of 
the Hamiltonian on the infinite strip can be obtained from the largest eigenvalue A G of the 
transfer matrix between different time slices as / = — In A G . To obtain A G , we used a Monte 
Carlo transfer-matrix methodic which has been shown to lead to accurate estimates of the 
largest eigenvalue even for models with continuous symmetry. The implementation of this 
method is similar to the case of the 2D frustrated classical XY model! and further details 
can be found in that work. 

The interfacial free energy for domain walls can be obtained from the differences between 
the free energies for the infinite strip with and without a wall. To obtain the critical expo- 
nents and critical coupling we employ the finite-size scaling AF(a, L) = A{L x l v ba) where 
A is a scaling function and 5a = a — a c . In a linear approximation for the argument of A, 
we have 

AF(a,L) = a + bL 1/u 5a (3) 

which can be used to determine the critical coupling a c and the exponent v independently!!!. 
The change from an increasing trend with L to a decreasing trend provides and estimate of 
a c . Once the critical coupling is known, the correlation function exponent rj can obtained 
from the universal amplitude a in Eq. (3) through a result from conformal invarianceii, 
a = 7T7]. To estimate the correlation length exponent v we first obtain the derivative S = 
dAF '/da near a c , then it can easily be seeing that a log-log plot of S vs L gives an estimate 
of \jv without requiring a precise determination of a c . In Fig. 2 we show this kind of plot 



for E x /E y = 3 from where we get the estimate v = 1.05(6). Of course, this is only valid 
in the linear approximation of Eq. (3). To ensure that higher-order terms can safely be 
neglected, the data for AF(a, L) must be obtained in a sufficiently small range near a c . We 
also checked, using a more general finite-size scaling analysisi, that allowing for higher-order 
terms gives results agreeing within the errors. The results for the critical coupling a c and 
critical exponents v and rj for two different values of the ratio E x /E y = 1 and 3 are indicated 
in Table I. 

For equal couplings E x = E y , the results for the critical exponents 77 and v differ sig- 
nificantly from pure 2D Ising exponents {y — 1 , 77 = 0.25). This result point to a single 
transition scenario. In fact, they are consistent with a point along the line of single tran- 
sitions in the XY-Ising modelEl. It is interesting to note that this result is also consistent 
with similar calculations for the 2D frustrated classical XY model!. In general, the critical 
behavior of a d dimensional quantum model is in the same universality class of the d + 1 
dimensional classical version. However, the ID FQXY model, apparently, is not the Hamil- 
tonian limit of the 2D classical model. Yet, their critical behavior appears to be in the same 
universality class. 

For the case of unequal couplings E x /E y = 3, the results indicated in Table I are in 
good agreement with pure 2D Ising values. From the relation between the ID FQXY model 
and the 2D XY-Ising model this then implies that the XY and Ising-like excitations have 
decouple in this region of parameters. To show that in fact one has two decoupled and, at 
the same time, separated transitions we show in Fig. 3 the results for the helicity modulus 
which measures the response of the system to an imposed twist. The helicity modulus is 
related to the free-energy differences AF between strips with and without and additional 
phase mismatch of 7r along the strip and is given by 7 = 2AF/-K 2 for large system sizes. 
If the model is decoupled then the transition should be in the universality class of the 2D 
classical XY model, where one knows that the transition is associated with a universal jump 
of 2/-7T in the helicity modulusS The critical coupling can be estimated as the value of a at 
which AF = 7T which gives a c = 1.29. This is to be compared with the critical coupling for 
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the destruction of chiral order in Table I, a c = 1.16. This clearly indicates the transitions 
are well separated and thus one expects they are decoupled. We have also performed less 
detailed calculations at other values of the ratio E x /E y from which we can estimate that the 
Ising and XY transitions merge into a single transition roughly at E x /E y ~ 2. Since, the 
superconductor to insulator transition is to be identified with the loss of phase coherence^ 
we reach the interesting result that in the ID FQXY, or alternatively, a Josephson-junction 
ladder, the universality class of the superconductor-insulator transition depends on the ratio 
between inter and intra-chain couplings. 
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FIGURES 

FIG. 1. Schematic representation of the one-dimensional frustrated quantum XY model with 
inter (E x ) and intra-chain (±E y ) coupling constants. The antiferromagnetic ordering of chiralities 
X p = ±1 is also indicated. 

FIG. 2. S = dAF(a,L)/da evaluated near the critical coupling a c . The slope of the straight 
line gives an estimate of \ jv. 

FIG. 3. Behavior of the interfacial free energy AF = L 2 Af for a system of size L = 12 
resulting from an imposed phase twist of ir. Vertical arrows indicate the locations of the Ising 
and XY transitions and the horizontal arrow the value AF = ir from where the XY transition is 
located. 



9 



TABLES 

TABLE I. Critical couplings (a c = (Ey/Ec) 1 / 2 ) and critical exponents (u, 77), obtained from 
finite-size scaling analysis of interfacial free energies for two values of the ratio between interchain 
and intrachain couplings (E x /E y ) 



E x /E y a c v 7] 

1 1.04(1) 0.81(4) 0.47(4) 

3 1.16(2) 1.05(6) 0.27(3) 
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